Astronomy 8l Astrophysics manuscript no. ms-pavan-igrJl 1 
July 15,2011 



©ESO2011 



IGRJ11014-6103: a newly discovered pulsar wind nebula? 

L. Pavan^ E. Bozzo^ G. Piihlhofer^, C. Ferrigno\ M. Balbo\ and R. Walter^ 



^ ISDC data center for astrophysics of the University of Geneva chemin d'Ecogia, 16 1290 Versoix Switzerland. 

^ Institut fiir Astronomic und Astrophysik, Kepler Center for Astro and Particle Physics, Eberhard-Karls-Universitat, Sand 1, 72076 



o 

(N 



X 

Oh' 

o 



> 

(N 

m 

00 
(N 

• 

O 



X 



Tubingen, Germany 
Received 2011 May 31; accepted 2011 July 09 



ABSTRACT 



Context. IGRJl 1014-6103 is one of the still unidentified hard X-ray INTEGRAL sources, reported for the first time in the 4th 
IBIS/ISGRI catalog. 

Aims. We investigated the nature of IGRJl 1014-6103 by carrying out a multi wavelength analysis of the available archival observa- 
tions performed in the direction of the source. 

Methods. We present first the results of the timing and spectral analysis of all the X-ray observations of IGRJl 1014-6103 carried 
out with ROSAT, ASCA, Einstein, Swift, and XMM-Newton, and then use them to search for possible counterparts to the source in the 
optical, infra-red, radio and y-ray domain. 

Results. Our analysis revealed that IGRJl 1014-6103 is comprised of three diff'erent X-ray emitting regions: a point-like source, an 
extended object and a cometary-like "tail" (~4 arcmin). A possible radio counterpart positionally coincident with the source was also 
identified. 

Conclusions. Based on these results, we suggest that the emission from IGRJl 1014-6103 is generated by a pulsar wind nebula 
produced by a high- velocity pulsar. IGR Jl 1014-6103 might be the first of these systems detected with INTEGRAL IBIS/ISGRI. 

Key words, stars: individual IGRJl 1014-6103 



1. Introduction 

The IBIS/I SGRI telescope (|Ubertini et all l2003h onboard 
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[aDSQQl 



INTEGRAL (J Winkler et al .1120031) has been operating since seven 
years from now, and the latest available catalog of sources 
detected ^vith th is telescope contains more that 700 objects 
(iBird et al.l l2QlQ). A large fraction of them belongs to the class 
of X-ray binaries (26%) and active galactic nuclei (35%); a rel- 
atively small fraction is represented by cataclysmic variables 
(5%), and another 4% comprises supernova remnants, globular 
clusters and soft y-ray repeaters. Abou t 30% of the IBI S/ISGRI 
sources are still unclassified (see e.g. lChaty et al.ll201Ql for a re- 
cent review). 

In this paper, we investigate the nature of the still uniden- 
tified INTEGRAL source IGR Jl 1014-6103. This source was re- 
ported for the first time in the 4th IBIS/ISGRI catalog (Bird et all 
|2Q1Q|) . and was detected in the IBIS/ISGRI mosaic at a sig- 
nificance level of 5.4 cr (20-100 keV). The best determined 
/A^reGTMLpositionof IGRJ11014-6103is at/?A=165.341 deg, 
Z)^c=-61.056deg (J2000), with an associated uncertainty of 
4.3 arcmin (90% c.L). The averaged source fluxes in the 20- 
40 keV and 40-100 keV energy band are 0.4+0.1 mCrab and 
0.6+0.2 mCrab, respectivel^Q. The ligh t curve of the sour ce, ob- 
tained from the online tool Heaveno Walter et aDl201Ql) . does 
not show evidence of variability. A possible counterpart to 
IGRJl 1014- 6103 in the soft X-ra y (0.3-10 keV) domain was 
identified by iMalizia et al.l (l201ll) . In the Swift /XRT field of 
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view (FOV) around IGRJl 1014-6103, they detected a source 
located at RA=165A4 deg, Dec=-61.022 deg (associated uncer- 
tainty ~6"). By using the refined X-ray position, these authors 
also suggested that IGR Jl 1014-6103 could be associated to the 
ROSAT source IRXH J110146.1-610121 and the serendipitous 
XMM-Newton source 2XMMJ1 10147.1-61012. 

In Sect. [2] we report on the analysis of all the available X-ray 
observations performed with Swift, XMM-Newton, ROSAT, and 
ASCA, in the direction of IGR Jl 1014-6103. We use the results 
of this analysis in Sect. [3] to search for possible counterparts to 
the source in the optical, infra-red and radio domains. A discus- 
sion on the nature of IGRJl 1014-6103 is presented in Sect.lH 

2. X-ray observations of IGR J11 01 4-61 03 

A summary of all X-ray observations analyzed in this section is 
provided in Table [T] 

2.1. Swift 

The FOV around IGRJl 1014-6103 was observed twice with 
Swift /XKT on 2011 March 10 and 12 for a total exposure time 
of 3. 1 ks and 2. 1 ks, respectively (observations ID. 00045395001 
and 00045395002). We analyzed the Swift /XRT (JGehrels et al.l 
2004) data collected in photon counting mode (PC) by using 
standard procedures (JBurrows et al.l 12005 ) and the latest cali- 
bration files available. Filtering and screening criteria were ap- 
plied by using ftools, and only event grades of 0-12 were 
considered. Exposure maps were created through the xrtex- 
POMAP task. No obvious sources were detected in the two ob- 
servations separately. We maximized the S/N by summing-up 
all the available data (eflfective exposure time 5.1 ks) and ex- 
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Fig. 1. t//?/7^r /^/i^: FOV around IGRJl 1014-6103 as observed by EPIC-MOSl on-board XMM-A^^w^^^ (0.5-12 keV). Upper right: 
Swift /XRT FOV around IGRJl 1014-6103 (1-9 keV). Bottom left: zoom of the EPIC-MOSl FOV. In this case we kept the colours 
and removed the coordinate grid and the error circles around the detected sources for clarity. Bottom right: zoom of the XRT FOV. 
In this figure we marked the source detected with Swift (uncertainty of 4.4" at 90% c.L), and sources N and S detected with XMM- 
Newton (uncertainty of 2" at 90% c.L). The contours determined from the XMM-Newton observations are also overplotted (in green), 
together with the error circle representing the INTEGRAL position of IGRJl 1014-6103 (in red) and the error circle of the Einstein 
source 2E 2383 (in blue, uncertainty of 39" at 90% c.L, see Sect 12.51) . The intensity scale, here and in the following images, is 
logarithmic and the grid is in galactic coordinates. 



Table 1. X-ray observation log of IGRJl 1014-6103. 



Instr 


OBSID 


START TIME 


Exp 






(MJD) 


(ks) 


Swift 


00045395002 


55632 


3.1 




00045395001 


55630 


2.1 


XMM-Newton 


0152570101 


52841 


60 




0111210201 


51753 


11 


ROSAT 


RH500445A01 


50627 


46 


ASCA 


51021010 


49733 


40 




51021000 


49412 


43 


Einstein 


2161 


44463 


11 



tracted a single image of the XRT FOV around IGRJl 10 14- 
6103. In this image, we detected a faint source (S/N=5.5 in the 1- 
9 keV energy band, see Fig.[T]) at the position 7?A= 165.444 deg, 
Z)^c=-61.023 deg (associated uncertainty 4.4" at 90% c.L; we 
used the xrtcentroid tool). This i s compatible with b eing the 
same source already identified by IMalizia et al.l (1201 lb . and is 



the only source detected inside the INTEGRAL error box. It 
is thus the most likely counterpart to IGRJl 1014-6103 in soft 
X-rays. Given the relatively short exposure time and the low 
count rate of the source, we could not extract a meaningful spec- 
trum. Instead, we estimated the source count rate (0.011+0.002 
counts/s, 1-9 keV) with sosta (available within the ftool xim- 
age), and used this result within webpi mms in order to de rive 
the corresponding X-ray flux (see also iBozzo et al.l l2009l) . In 
the conversion between count rate and flux, we adopted an 
absorbed (A^h=0.7x10^^ cm"^) power law (r=1.6) model (see 
Sect. 12.21) . The inferred X-ray flux in the 2-10 keV energy band 
was (6.3±1.2)xl0"^^ erg cm"^ s"^ (not corrected for absorp- 
tion). Even though a single point source is detected with the 
XRTCENTROID tool. Fig. [T]reveals a marginal evidence for difl'use 
emission around the source and the presence of an elongated 
structure (~4 arcmin) extending in the north-west direction from 
the point source. We discuss further these evidences in Sect 12.21 
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Fig. 2. Left: close view of the EPIC-MOS 1 FOV around sources 
N and S (see Sect.O; Right: simulated PSF of the EPIC-MOSl 
camera at energies in the (0.2-12) keV energy range by assuming 
the same off-set position as that of sources N and S. The two 
panels have the same angular scale. The simulated PSF does not 
show the double peak profile seen in the real data. We also note 
that the distortion due to the off'-set position of the two sources 
would occur in a diff'erent direction with respect to the separation 
between the two sources. 



2.2. XMM-Newton 

XMM-Newton observed the region around IGRJl 1014-6103 
on 28 July 2000 (obs. ID 0111210201, total exposure time 
-11 ksec) and on 21 July 2003 (obs. ID 0152570101, total expo- 
sure time ~ 60 ksec). In both cases, the observations were cen- 
tered on the nearby supernova remnant (SNR) MSH 1 1-61 A, and 
thus IGRJl 1014-6103 was off'-set with respect to the instrument 
aim-point by about 15 arcmin. During obs. ID 0152570101, the 
EPIC-pn and EPIC-MOS cameras were operated in small win- 
dow and full frame mode, respectively. IGRJl 1014-6103 was 
only within the EPIC-MOS FOV. During obs. ID 0111210201 
all the EPIC cameras were operated in full frame mode, and 
IGRJl 1014-6103 was observed with both the EPIC-pn and 
EPIC-MOS. In all cases, only part of the INTEGRAL error circle 
around IGRJl 1014-6103 was included in the FOV of the EPIC 
cameras (see Fig.[T]). 

We processed the XMM-Newton observation data files with 
the two pipelines epproc and emproc in order to produce EPIC-pn 
and EPIC-MOS cleaned event files, respectively (SAS v. 10.0.1). 
To identify the high background time intervals, we followed the 
SAS online analysis threacfl and extracted the lightcurves from 
the full FOV of the EPIC cameras in the 10-12 keV energy band. 
We discarded from further analysis the time intervals in which 
the 10-12 keV FOV count rate was larger than 0.4 cts/s. The 
resulting eff'ective exposure times were 47.3 ks and 8.0 ks (5.7 
ks) for the EPIC-MOSl (EPIC-pn) cameras in the observation 
ID 0152570101 and 01 11210201, respectively. All cleaned event 
files were barycentered with the barycen tool. 

A visual inspection of the EPIC-MOS image extracted from 
the XMM-Newton observation with the longer exposure time 
(ID 0152570101, see Fig. [T]) reveals that the source detected 
with Swift jXKV is indeed composed of three diff'erent emit- 
ting regions: two sources separated by ~ 22 arcsec (hereafter 
sources N and S) and a dimmer elongated structure (~4 ar- 
cmin, hereafter, the "tail"). In the EPIC images extracted from 
obs. ID 0111210201 only the two brighter sources are clearly 
visible. In this observation, the tail is only marginally detectable 
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Fig. 3. Zoom of the EPIC-MOSl FOV around IGRJl 1014-6103 
during observation ID. 0152570101 (0.2-12 keV). The regions 
used for the extraction of the spectra of sources N and S are 
indicated with red small circles; the region used for the spectral 
analysis of the tail with a red ellipse. We also show the extraction 
region used to perform the spectral analysis of the entire emitting 
region (dashed ellipse; see Sect. 12.2b . 



^ See also http://xmm. esac.esa.int/sas/current/documentation/threads/ 1 
PN_spectrum_thread. shtml. 



due to the shorter exposure time. In order to check that sources N 
and S are not just produced by the distortion of the EPIC-MOS 
point spread function (PSF) at large off'-set angles, we simulated 
with the psfgen task the 2D-PSF of the instrument at the location 
of these two sources. We used a circular region for the simula- 
tion -180 arcsec wide centered on source N and set the level pa- 
rameter to "ELLBETA". The comparison between the simulated 
PSF profiles and the measured intensity profiles obtained from 
the data is shown in Fig. O The simulated PSF did not present 
a double peak in the intensity profile in a direction comparable 
with that observed from the XMM-Newton data. In particular, we 
note that the distortion of the PSF caused by the off'-axis posi- 
tion of the source on the EPIC-MOS detectors occurs preferably 
along the azimuthal direction, whereas sources N and S are sep- 
arated in the opposite direction. 

Sources N and S were also already identified in 
the sec ond XMM-Newton catalog of serendipitous X-ray 
sources (IWatson et al.l l2009t) as 2XMMJ110147.1-610124 and 
2XMMJ1 10145.0-610140, respectively. The best determined 
positions are (/?A= 165.4465 deg, Z)^c=-61.0234deg) for the 
first, and (7?A= 165.4377 deg, Z)^c=-61.0279deg) for the lat- 
ter (the nominal positional accura cy of the EPIC cameras 
is ~2", see e.g. iPavan et aP l201ll and references therein). 
2XMMJ110147.1-610124is reported as extended, and is char- 
acterized by a radial elongation of 8.1 arcsec. 

In Fig. [T] the image extracted from Swift jXKV is compared 
with that obtained from the deeper XMM-Newton observation. 
The contours around the extended emission from XMM-Newton 
coincide well with the emission in XRT. This confirms the 
marginal detection of the extended emission noticed before in 
the Swift IXKV data (see Sect. 12. lb , and suggests that the single 
point source detected with this telescope is indeed a blend of the 
two brighter objects revealed by XMM-Newton. 

We also performed a spectral analysis of the X-ray emis- 
sion from the two brighter sources and the tail extracting 
the corresponding spectra from the EPIC-MOSl data in obs. 
ID 0152570101 (EPIC-M0S2 data could not be used for the 
spectral analysis as in this camera the sources were located 
close to the gap between two diff'erent CCDs). The extraction 
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Fig. 4. Contour plot of the absorption column density (Nu) and 
power law photon index (F) derived from the EPIC-MOS 1 data 
for sources N and S in observation ID. 0152570101. The con- 
tours correspond to 68%, 90% and 99% confidence levels. 



regions (shown in Fig. [S]) were chosen in order to maximize 
the S/N and to minimize the reciprocal contamination. To take 
into account the background contribution, we selected a circular 
source-free region of radius 2 arcmin located in the same CCD as 
IGRJl 1014-6103 and at approximately the same distance from 
the central SNR. We also checked that the results presented in 
this paper would not be aff'ected by diff'erent (reasonable) choices 
of the background region. Given the relatively low number of 
counts, all the spectra were grouped in order to have at least 5 
photons for each energy bin and fit with an absorbed power law 
model using C-statistics (Cash 1979, ; we used Xspec v. 12). We 
checked that similar results (to within the uncertainties) would 
have been obtained by grouping the spectra with 25 photons per 
energy bin and using ;^^ statistics. The results of the fits are sum- 
marized in Table O Here we also report the best-fit parameters 
obtained from the spectral analysis of the entire region (dashed 
ellipse in Fig. [3]). We found marginal evidence for a diff'erence 
between the power law photon index of sources N and S, and 
show in Fig. m the most relevant contour plots derived from the 
spectral fits. The spectrum extracted from the tail provided only 
poor constraints on the properties of its X-ray emission. A sim- 
ilar analysis carried out on the EPIC-MOS and EPIC-pn data of 
the observation ID. 0111210201 did not reveal any significant 
variation of the estimated spectral parameters and fluxes from 
the three regions. 

In order to compare the fluxes estimated with XMM-Newton 
and those reported for IGRJl 1014-6103 in the 4th IBIS/ISGRI 
catalog (see Sect. [B, we performed a simultaneous fit of the 
MOS 1 (ID. 0152570101) spectrum from the entire source 
(dashed ellipse in Fig. O and the IBIS/ISGRI spectrum of 
IGRJl 1014-6103 obtained from the Heavens web-tool. The two 
spectra could be reasonably well fit using an absorbed power law 
model, with spectral parameters compatible with those deter- 
mined previously using only the XMM-Newton data (best-fit val- 
ues reported in the last line of table[2l). We obtained A^h=0.6±0.2, 
F= 1.5+0.2 OtVd.o.f. =0.8/74) and C=1.7:^;^. Here, C is the 
normalization constant introduced to take into account both 
the intercalibration between the EPIC-MOS 1 and IBIS/ISGRI, 
and any possible flux change of compact emission regions be- 
tween the difl'erent observations. We considered EPIC-MOS 1 
as reference detector, therefore the constant is applied to the 



Table 2. Results of the spectral fit parameters obtained from 
XMM-Newton observation ID. 0152570101. AU spectra were 
fit using an absorbed power law model (photon index F). The 
absorption column densities (A^h) and fluxes are in units of 
10^^ cm"^ and 10"^^ erg cm"^ s"\ respectively. We also report 
the estimate of the X-ray flux in the 0.2-2.4 keV energy band 
(^0.2-2.4 kev), in order to have an easier comparison with the re- 
sults obtained from ROSAT/HRl observations (see Sect. 12.3b . 
Uncertainties are at 90% c.l. on the spectral parameters and 68% 
c.l. on the fluxes. 





A^H 


F 


^2-10 keV 


Fo.2-2A keV 


C-stat/d.o.f. 


N 


0.7+0.2 


2.0+0.3 


3 7+0.6 


1-2!^:? 


68.2/68 


S 


1.0+0.4 


1.3+0.4 


"•^-2.6 


o-ni 


73.3/62 


Tail 


1 3+^-^ 


9 9+0.7 
^•^-0.8 


1 Q+0.3 
^•^-1.0 


^"^-0.3 


31.1/21 


Total 


0.7 + 0.1 


1.6 + 0.15 


^7+1.2 


o Q+0.3 
•^•^-0.5 


68.60/85 



ISGRI spectrum. We noticed that fixing C=l would not signif- 
icantly change the results of the fit. The estimated 2-10 keV, 
20-40 keV, and 40-80 keV X-ray fluxes of the sources were 
(1.9 + 0.2)xl0-i2 ergcm-2s-^ (3 A^_^^j)xlO-^^ erg cm'^ s-^ 
and (4.8^Q^)xlO"^^ erg cm"^ s"^ respe ctively. Th ese are fully 
compatible with those reported by Bird e t al.l (1201 0), i.e. (3.0 + 
0.8)xlO-i2 and (5.6 + 2)xl0-i2 erg cm'^ s'^ in the 20-40 keV 
and 40-100 keV band respectively, see Sect. [T] The unfolded 
MOSl+ISGRI spectrum of IGRJl 1014-6103 is shown in Fig.[5l 
A simultaneous MOSl+ISGRI spectral analysis was carried out 
also using the EPIC spectra of source S alone (without includ- 
ing the tail and source N). This source indeed seems to have an 
harder spectrum and higher estimated flux in the 2-10 keV band 
with respect to source N and the tail. No significant diff'erences 
in the three spectral parameters have been measured, except for 
a marginal hint of an increase of the intercalibration constant 
needed in the fit. 

We also searched in the EPIC data for indications of vari- 
ability for sources N and S on time scales of seconds to hours, 
but did not spot any relevant feature. We paid particular atten- 
tion in searching for coherent pulsations. We used the PN data 
for the short observation (ID. 0111210201) and decided to ana- 
lyze sources N and S jointly, to maximize the number of source 
photons. Cutting out source S for an analysis of source N (and 
vice versa) would result in a too high loss of source photons for a 
detailed timing analysis, respectively. The total number of pho- 
tons obtained for both sources together was 693. By adopting 
an epoch-folding technique with 16 phase bins, we determined 
an upper limit on the pulsed fraction for a sinusoidal signal of 
56% (at 90% c.l.) in the 0.002 to 6.8 Hz fr equency range (we 
followed the approach of iLeahy et al.lll983h . In the long obser- 
vation (ID. 0152570101) we could use only the data of MOSl 
(1275 photons) and M0S2 (2087 photons), and explored there- 
fore the frequency range 0.001-0.2 Hz: the upper limits on the 
pulsed fraction are 44% and 33% (at 90% c.l.) respectively. 

2.3. ROSAT 



We retrieved all available ROSAT (iTruemped Il982h data in 
the direction of IGRJl 1014-6103 from the Heasarc archive. 
The FOV around IGR J l 1014 -6103 was observed with the 
HRI (Pfeff'ermanneLalJ Il987h on-board ROSAT on several 
occasions (observations ID. RH500445A01, RH500340N00, 
RH500445N00, RH900619A01, and RH900619N00). Among 
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Fig. 5. Unfolded XMM-Newton and IBIS/ISGRI spectrum of 
IGRJl 1014-6103. We fitted both spectra with a single absorbed 
power law model (solid line). In the bottom panel we show the 
residuals from this fit. 



these, observation ID. RH500445A01 was characterized by the 
longest available exposure time (46 ks) and was performed on 
1997 June 28. We show in Fig.[6]the HRI image (5" resolution, 
0.2-2.4 keV). Within the INTEGRAL error circle of IGRJl 1014- 
6103 a single point-like source is detected. This source was al- 
ready classifi ed in the HRI catalog as IRXHJl 10146.1-610121 
("Voges etaP 119991) . The best determined source position is at 
7?A= 165.442 deg and Dec=-61.022 de^. The nominal positiona l 
accuracy of the HRI is of order 10" (iPredehl & Prietol l200lh 
and thus the ROSAT position is consistent with that of the 
source revealed with Swift and with source N detected with 
XMM-Newton. Given the limited spatial resolution, sensitiv- 
ity and energy coverage of the HRI with respect to the EPIC 
cameras, it is not surprising that the closeby source S and the 
tail are not detected in this case (source S is also fainter than 
source N in the ROSAT energy band, see Table [2l). We sug- 
gest that, as for the case of Swift /XKT, the ROSAT source is in- 
deed a blend of sources N and S. The estimated count rate of 
IRXHJl 10146.1-610121 in the HRI energy band (0.2-2.4 keV) 
is 0.0038+0.0004 cts/s. This would convert into an X-ray flux of 
(9±l)xl0"^^ erg cm"^ s"^ (2-10 keV) if an absorbed power law 
spectral model with A^h=0.7x10^^ cm"^ and r=1.6 is used (see 
Sect. 12.21) . The flux of the ROSAT source is thus qualitatively in 
agreement with that of sources N and S estimated from XMM- 
Newton data. 

We noted that the FOV around IGRJl 1014-6103 was ob- 
served also with the PSPC on-board ROSAT on two occasions, 
on 1993 July 23 (observation ID. RP500307N00) and on 1994 
July 5 (observation ID. RP500307A01). However, in these two 
cases the relative low exposure time (~2 ks) and the ofl'-set po- 
sition of the source with respect to the instrument aim-point, did 
not result in a clear detection of the source. We thus do not dis- 
cuss these data in more detail. 
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Fig. 6. Top: ROSAT /URl FOV around IGRJl 1014-6103 (obser- 
vation ID. RH500445A01). The image is in the 0.2-2.4 keV en- 
ergy band with a resolution of 5" (see Sect. 12.31) . We overplot on 
this image the XMM-Newton contour levels (in blue) determined 
from the EPIC-MOSl data in observation ID. 0152570101 (see 
also Fig. [B. The black circle represents the INTEGRAL posi- 
tion of IGRJl 1014-6103. The bright extended emission in the 
center of the FOV is SNR MSH 11-61 A. Bottom: A^CA /GIS2 
FOV around IGRJl 1014-6103 (observation ID. 51021000). The 
image is in the 0.5-10 keV energy band. Extended emission is 
clearly visible in the bottom right side of the detector, below 
SNR MSH 1 1-61 A. We overplot on this image the XMM-Newton 
contours (in white), the INTEGRAL error circle of IGRJl 1014- 
6103 (red circle), and the extraction region adopted for the GIS 
spectral analysis (yellow ellipse). 



2.4. ASCA 

In the Heasarc archive we fo und two observatio ns carried 
out wi th the GIS telescope (" Ohashi et al.l [1996') on-board 
ASCA ( Tanaka et al.l 1994) which included the region around 
IGRJl 1014-6103 in their FOV. The first of these observations, 
ID. 51021000, was performed on 1994 March 1 with a total ex- 



posure time of 43.4 ks. The second observation, ID. 51021010, 
lasted for 40.5 ks and was carried out about 10 months later, 
on 1995 January 16. Both these observations were aimed at the 
SNR MSH 1 1-61 A and did not cover entirely the INTEGRAL er- 
ror circle around the position of IGRJl 1014-6103. However, the 
extended emitting region already identified with XMM-Newton 
is clearly detected in the GIS images (see Fig. [6]). In comparison 
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with the HRI camera on-board ROSAT, the GIS has a broader 
energy coverage (0.5-10 keV) and a larger sensitivity, but a sig- 
nificantly worse spatial resoluti on (the GIS P SF is ~ 1 arcmin and 
strongly dependent on energy. Tanaka et aD il994). It is thus not 
surprising that the GIS is not able to distinguish the two sources 
(N and S) detected with XMM-Newton and ROSAT, but can re- 
veal the presence of the diff'use emission around them. 

We analyzed GIS 2 and GIS 3 data following standard pro- 
cedures and the latest available calibration files (see e.g. 
iPavan et aP l201l l and references therein). Source spectra in 
the two observations were extracted using the elliptical region 
shown in Fig. [6l Given the presence of the nearby SNR MSH 
11 -61 A, we extracted the background spectra from diff'erent 
regions located at the same radial distances of IGRJ11014- 
6103 from MSH 11-61 A. Even though the S/N of the data 
was relatively poor, the GIS spectra could be reasonably well 
fit using an absorbed power law model. From the GIS2 data 
in observation 51021000, we measured r= 1.6+0. 3 and esti- 
mated a 2-10 keV X-ray flux of (2.1+^-2)xl0-^2 erg cm'^ s"^ 
(X^^ J d.o.f =0.1/21). These values are compatible (to within 
the errors) with those inferred from the XMM-Newton data. 
From observation 51021010, we obtained r= 1.2+0. 5 and 
/^2-iokeV=(2.1!?i)xlO-i2 ergcm-2s-i Orred/d.o.f.=0.9/16). In 
both cases we fixed the absorption column density to the best 
value measured from XMM-Newton. Compatible results (to 
within the uncertainties) were obtained from the GIS 3 data. 
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Fig. 7. Search for optical and infra-red counterparts to the 
sources N and S. Left: STCI-DSS image. We overplot the po- 
sition of sources N and S (radius equal to the position accu- 
racy, 2") as determined from the EPIC-MOSl observation, and 
the most likely optical counterpart from the USNO-A2.0 and 
USNO-Bl.O catalogs (green circles, with radius 0.6"). Right: 
2MASS image. The position of sources N and S is indicated as 
above. We marked the location of the potential infra-red coun- 
terparts retrieved from the 2MASS catalog (green circles, with 
radius 0.6"). 



2.5. Einstein 

The region around IGRJl 1014-6103 wa s observed also with 
the Einstein/IFC detector (0.2-4.5 keV, iGiacconi et sl\ Il979h 
on 1980 August 12 with an exposure time ~ llks. A vi- 
sual inspection of the Einstein/IFC image revealed a source 
(2E 2383) spatially coincident with sources N and S and a 
marginal evidence for an extended tail similar to that observed 
in XMM-Newton. The best position of 2E 2383 is reported 
in the 2E catalog at 7?A = 11^ : 01"^ : 47.5\ Dec= - 6r:0y :2r 
(related uncertainty 39'', see FigUl iHarris et al.l Il994l) . The 
source count rate was 0.011+0.002 cts/s, which correspond to 
(9 + 2) X 10"^^ erg cm"^ s"^ (assuming a column density Nh = 
0.7 X 10^^ cm"^ and a power law spectrum with slope F = 1.6, 
see Sect. 12.2b . This flux is compatible with the one measured 
with XMM-Newton. 



3. Counterparts to IGRJ11014-6103 

In this section we take advantage of the improved positions for 
the three emitting regions comprised in the INTEGRAL error 
circle around IGRJl 1014-6103 to search for the correspond- 
ing counterparts in the optical, infra-red, and radio domain. In 
Sect. 13.31 we also report on a possible y-ray detection. 



3.1. Optical/1 R data 

We used publicly available optical and near infra-red images 
from STSCI-DSS, and 2MASS surveys. Potential counterparts 
to sources N and S are shown in Fig. [71 For source N, we 
did not find any obvious counterpart in the optical domain. 
The closest catalogued USNO-A2.0 sources are located more 
than ~7" away from the best determined EPIC-MOSl position 
of the source. The nearest classified object in the USNO-Bl.O 
catalog is J0289-0254424 (Rl = 10.89, B2=15.97), located at 
/?A= 165.4471 deg, Z)^c=-6 1.0221 deg. In the infra-red domain. 



the closest catalogued source is the 2MASS object Jl 1014730- 
6101213 located at 7?A= 165.447 122 deg, D^c— 61.022602 deg. 
This object is reported in the catalog with J= 15.95, H= 15.54, 
K= 15.46. For source S, the closest catalogued source in 
the optical domain (USNO-A2.0 catalog) is J0225- 10220547. 
This is located at 7?A= 165.438870 deg, Z)^c=-6 1.027375 deg 
and is characterized by B=16.5 and R=13.0. In the infra- 
red domain, the closest object from the 2MASS catalog to 
the position of source S is Jl 1014532-6101383, located at 
/?A= 165.438860 deg, Z)^c=-6 1.027306 deg and characterized 
by J=11.16, H=10.19, K=9.89. 

3.2. Radio data 

To search for radio counterparts to IGRJl 1014-6103, we used 
data from the MGPS-2 archive. The MGPS-2 is a high- 
resolution and large-scale survey of the galactic plane car- 
ried out with the Molonglo Observatory Synthesis Te l escop e 
(MOST) at a frequency of 843 MHz (iMurphv et al.l l2007h . 
The MGPS-2 image covering the field around IGRJ11014- 
6103 is shown in Fig. [8l We noticed from this figure the 
presence of a relatively bright radio source, Jl 10149-610104, 
positionally consistent with sources N and S. The radio 
source is reported in the MGPS catalog as "compact", be- 
ing its dimensions comparable to the convoluted beam of the 
survey. The best determined position of Jl 10149-610104 is 
7?A=(165.4559+0.0017)deg, Z)^c=(-61.0178+0.0017)deg. The 
associated flux at 843 MHz is 24.2+4.8 mJy (corresponding to 
(2.04+0.40) xlO-^6 erg cm'^ s"^). 



3.3. y-ray data 

We also noticed the presence of an EGRET source, 3EG Jl 102- 
6103, which might be related with IGRJl 1014-6103. In the lat- 
est available EGRET catalog, the best determined position of 
3EG J1102-6103 is at /?A- 165. 60 deg, Z)^c=-6 1.05 deg, with 
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Fig. 8. Left: radio image (and contours) at 843 MHz around 
the position of IGRJl 1014-6103 retrieved from the MGPS- 
2 archive. We mark in red the position of the source 
MGPS2J1 10149-610104. Right: EPIC-MOSl image from ob- 
servation ID. 0152570101. We overplot on the image the radio 
contours derived from the MGPS-2 data in green and marked in 
red the position of the radio source Jl 10149-610104 and the two 
XMM-Newton sources N and S. 



an associated uncertainty of 0.6 deg. The error circle around 
3EG Jl 102-6103 would thus be compatible with the position of 
IGRJl 1014-6103, as determined before. At present, however, 
we are unable to confirm this association due to the presence 
of other possible counterparts located within the EGRET error 
circle of 3EG Jl 102-6103. Among these, there are the SNR 
MSH 11 -61 A and the young energetic pulsar PSR Jl 105-6107 
(lKaspietal.1119971) . 

In the Fermi/LAI catalog we did not find any source posi- 
tionally coincident with IGRJl 1014-6103. We did not attempt 
an estimate of an upper limit on the source flux from the Fermi 
data. The latter calculation would indeed depend crucially on 
the properties of the still relatively poorly known high energy 
Galactic background in such a highly crowded region. 

4. Discussion and conclusions 

In this paper we analyzed all publicly available X-ray ob- 
servations carried out in direction of the unidentified source 
IGRJl 1014-6103 with Swift, XMM-Newton, ROSAT, ASCA and 
Einstein. The entire INTEGRAL error circle around IGR Jl 1014- 
6103 has been covered in the soft (1-10 keV) X-ray energy band 
with Swift and ROSAT. The analysis of these data led to the de- 
tection of a single X- ray source inside the error circle given by 
INTEGRAL (see also iMalizia et al. l '2011*). However, the deep- 
est observations of the region around IGRJl 1014-6103, per- 
formed with the EPIC cameras on-boardXMM-A^^w^6>^, revealed 
that the source detected with Swift/XKT and ROSAT is likely a 
blend of two distinct objects: a point-like source and an 8.1 arc- 
sec extended object. The same observations revealed also the 
presence of a cometary-like tail (~4 arcmin) extending from the 
two sources, and comprised within the INTEGRAL error circle 
around IGRJl 1014-6103. The X-ray spectra of the three objects 
displayed only marginally significant diff'erences (see Table |2] 
and Fig. [4]). 

The ROSAT /URl and ASCA/GIS observations of the same re- 
gion were unable to disentangle the three emitting components. 
This can be ascribed to the narrower energy range coverage of 



the former instrument and the lower spatial resolution of the tele- 
scopes on board ASCA and ROSAT with respect to the EPIC 
cameras (see fluxes in the 0.2-2 A keV band in Table [2] and 
Sect. 12.31 12.41) . In all cases, we showed that the X-ray emitting 
properties of the entire region comprised in the INTEGRAL er- 
ror circle of IGRJl 1014-6103 were in agreement among all the 
available instruments. In Sect. 12.21 we also showed that the joint 
XMM-Newton-\-INTEGRAL spectrum could be reasonably well 
described using a single absorbed power law model, with a nor- 
malization constant between the two instruments ^ 1 . All these 
results suggest that IGRJl 1014-6103 is a persistent X-ray emit- 
ter that displayed in the past -31 years only marginal variations 
(if any) in the 0.3-100 keV energy band. 

The improved X-ray positions of sources N and S per- 
mitted also to search for their possible optical and infra-red 
counterparts. The closest catalogued 2MASS, USNO-A2.0, and 
USNO Bl.O objects to sources N and S lie just outside the 90% 
c.l. error circle determined with XMM-Newton; further observa- 
tions with the ACIS telescope on-board Chandra are therefore 
needed to reduce the likelihood of a chance coincidence. We 
found a possible radio counterpart close to sources N and S. 
The radio contours derived from the publicly available MGPS-2 
data match reasonably well with those of the X-ray emission 
detected with XMM-Newton around sources N and S. The 
presence of an EGRET source, 3EG Jl 102-6103, positionally 
coincident with IGR Jl 1014-6103 was also reported. 

Being IGRJl 1014-6103 located only 11 arcmin apart from 
the nearby SNR MSH 11 -61 A, we first considered the pos- 
sibility that the X-ray emission coming from the unidentified 
INTEGRAL source was related to the SNR. However in none 
of the investigated wavelengths (X-ray, radio, optical, NIR) we 
could clearly identify the presence of emitting structures located 
between MSH 11-61 A and IGR Jl 1014-6103 which could have 
favored such an association (see also Rosado et al. 1996). 

As an alternative and perhaps more viable hypothesis, 
we considered that the morphology and emission properties 
of the 4 arcmin extended X-ray tail in IGRJl 1014-6103 
might resemble the elongated feature s observed in the ca se of 
PSR B2224+65 (" the Guitar" nebula. lHui& Becked l2007h and 



PSR J0357+3205 (JDe Luca et aDl201 ih . In both cases a pulsar 
(PSR) has been firmly detected at one end of the elongated struc- 
tures, and the cometary-like tails have been tentatively explained 
under the "bow-shock pulsar wind nebulae" (bsPWN) scenario. 
These elongated cometary-like objects form when a pulsar es- 
capes from the associated supernova remnant and moves with 
a very high velocity through the interstellar medium (typically 
hundreds of km/s, i.e. orders of magnitude higher than the ve- 
locities of the col d and warm com ponents of the interstellar 
medium, see e.g. Roberts et al.l 12005 . and references therein). 
Close to the pulsar, the relativistic wind of the compact object 
usually gives rise to a sub-luminous cavity in X-rays surrounded 
by a termination shock. In this region, particles are thermal- 
ized and accelerated, thus producing a conspicuous X-ray and 
radio emission (through synchrotron processes due to the inter- 
action with the local magnetic field). As the pulsar is moving, 
a bow shock is formed in front of it and the flow of material 
advects the emitting particles back along the direction of mo- 
tion, lea ding to the formation of a n elongated cometary structure 
(see e.g. Gaensler & Slanei 120061 for a detailed discussion). In 
a few cases, deep X-ray observations with the ACIS telescope 
on-board Chandra permitted to clearly disentangle the difl'erent 
contributions to the total X-ray and radio emission produced by 
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the structures around a bsPWN (se e e.g. the case of the PWN 
G359.23-0.82;'GaenslereLaDl2QQ£). 

If the structures detected in IGRJl 1014-6103 are indeed 
originating from the pulsar wind, the X-ray tail can be inter- 
preted as being due to relic electrons spread along the direc- 
tion of motion of the PSR. Source N, extended in X-rays and 
closer to the peak of the radio emission, would represent, ac- 
cording to this interpretation, the compact PWN (with possi- 
ble contribution coming also from the bow- shock). The detec- 
tion of a radio source located close to source N would sup- 
port this scenario, as radio synchrotron emission is expected 
from particles accelerated a t the termination s hock; see th e cases 
of, e.g. PSRs B1853+01 (iFrail et al.1 [T995 iPetre et aP [2002) 
B1957+20_(StaDDers eraD^Sr^757-24 (iFrail & Kulkarnil 
ll991l:lKaspietal.n200Th . 

While the association between the extended source N ad 
the "tail" seems morphologically plausible, the association of 
source S to this region is less straightforward. 

In case source S is not physically related to the X-ray tail, 
it might just by chance be located along the line of sight to the 
extended structure. The pulsar responsible for the formation of 
the compact and relic PWNe, might, in this case, be unresolved 
and comprised in the emission detected from source N. 

In case, instead, that source S, the point-like source, is phys- 
ically connected to the PWN, it would be reasonable to associate 
it to the PSR generating the X-ray structures. Under this scenario 
the misalignment between the direction of the 4 arcmin tail and 
the axis between sources N and S would imply a significant tilt 
between the proper motion of the PSR and the position of the 
compact PWN. A similar misalignment, however, was already 
observed in e.g. PSR B2224+65, and ascribed to the strong in- 
teraction betwee n the fast-moving p ulsar and its dense surround- 
ing environment (iHui & Becke i^ 2007). Indeed, as for the case of 
PSR B2224+65, also the environment around MSH 1 1-61 A and 
IGRJl 1014-6103 was found to be a complex and dense region, 
characterized by the prese nce of several backg round and fore- 
ground molecular clouds (iFilipovic et al .1120051 and references 
therein). 

Following the joint XMM+ISGRI spectral analysis (see 
Sect. 12.21) , we conservatively assume that both sources N and S 
contribute to the ISGRI emission (the tail component has indeed 
a much lower flux). Spectral properties (hardness and flux sta- 
bility) similar to those observed in the present case (see Fig. [5]) 
have already been measured in, e.g. AX Jl 838-0655. This 
source is a young (r^ = PI2P = 2 3 kyr) pulsar, surrounded by a 
PWN (iGotthelf & Halpernll2008h . The best-fit po wer law slope 
in the 2-100 keV was measured at F = 1 .5 + / - 0.2 (iMalizia et al.l 
l2005h . We note however that also the highly-magnetized neutron 
stars C'magnetars", for a review see e.g. Mereghetti 2008) share 
similar spectral properties. SGR 1806-20, for example, showed 
a power law s lope F = 1.6-2.0 in the range 1-100 keV 
(lEsposito et al.l 12007). As the 1-10 keV spectra of magnetars 
are usually (phenomenologically) modelle d with a combined 
black body+power law (BB+PL) model (M ereghetti2008 ). the 
lack of the BB component in the XMM-Newton spectrum might 
be an issue for this interpretation. We therefore estimated for 
source S (i.e. the component showing the higher 2-10 keV flux) 
the contribution of an eventual black body spectral component. 
We used a temperature kT = 0.5, in agreement wit h the one 
measured for SGR 1806-20 (lEsposito et al.1 l2007h and for 
magnetars in general. Due to the relative low S/N, we could only 
obtain a poorly constrained upper limit on the ratio between 
the fluxes of the BB and PL components /^bb/^pl < 0.1. This 
would be compatible with the value measured in the case of 



SGR 1806-20 (Fbb/Fpl = 0.03, lEsposito et all l2007h . The 
present data therefore do not argue against a magnetar hypoth- 
esis. We note also that a PWN associated to a magnetar has 
already been observed , e.g. in the case of AXP 1E1547. 0-5408 
dVink & Bamball2009l) . 

Further observations with an X-ray telescope characterized 
by a finer spatial resolution (i.e. the ACIS on-board Chandra) 
and increased spectral and timing statistics are required to 
solve the issues above and firmly establish the real nature of 
IGRJl 1014-6103. In case these observations will confirm that 
IGRJl 1014-6103 is a newly discovered PWN generated by 
a high-velocity PSR, this would be the first detection with 
INTEGRAL of one of these systems (to the best of our knowl- 
edge). 
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